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Supernovae provide a unique opportunity for large scale parity violation because they are domi- 
nated by neutrinos. We calculate the parity violating asymmetry A of neutrino emission in a strong 
magnetic field. We assume the neutrinos elastically scatter many times from slightly polarized neu- 
trons. Because of the multiple interactions, A grows with the optical thickness of the proto-neutron 
star and may be much larger than previous estimates. As a result, the neutron star could recoil at 
' a significant velocity. 
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Weakly interacting neutrinos dominate core collapse supernovae. This provides a unique opportunity for large 
, scale or macroscopic parity violation [1,2]. Therefore, it is of fundamental importance to study parity violation in a 
supernova. 

We focus on an asymmetry induced in the explosion because of parity violation in a strong magnetic field. This 
' asymmetry could lead to a recoil of the newly formed neutron star. Indeed, neutron stars have large velocities ~ 500 
km/s [3]. An asymmetry of order one percent could produce these velocities [4] (see below). However, this asymmetry 
need not arise from parity violation (sec for example [5]). 

Strong magnetic fields are present in pulsars. Indeed, external dipole fields of 10^^ to 10^"^ G are inferred in many 
' cases [6]. In this paper, we estimate the magnitude of parity violating effects from known neutrino interactions in 
CsJ known fields of 10^^ to lO^'^ G. Others have speculated on parity violating effects in much stronger fields [1,2,7-11] 
and with new neutrino interactions [12]. We find that repeated interactions as neutrinos diffuse through an optically 
thick medium may greatly enhance the asymmetry. 

There is some controversy on observational correlations between recoil velocities and magnetic field strength. For 
, example, Birkel and Toldra [13] argue that for rapidly spinning pulsars there is no correlation between the recoil 
velocity and the projection of B on the spin axis. We think their analyses may be over simplified because they do 
not consider possible effects of rapid rotation on the dynamics of the collapse and on the asymmetry. Unfortunately 
observational tests involve incomplete information. The strength of the external dipole field is inferred. However, 
little is known about non-dipole fields. We find that the most important variable may be the volume of the core 
occupied by the strong field (see below) rather than simply its strength. Thus we keep an open mind with respect to 
^ , present observations. 

■ Independent of observation, it is important to estimate parity violating effects. We see three possibilities: parity 
• I— I , violating effects could be small and thus irrelevant, they could be large and observed or they could be potentially 
' large and not observed. If they are not observed, it may still be possible to set useful limits on the magnetic field 
configuration or on new weak interactions. In any case, we need accurate theoretical estimates. 

Previous estimates of parity violation may be incomplete because they ignore possible enhancements from repeated 
interactions. Neutrino transport involves diffusion with neutrinos undergoing many parity violating interactions before 
they escape. We estimate these cumulative effects below. 

Much previous work focused on electrons [1,11]. It is natural to think that neutrino electron scattering will dominate 
the asymmetry because the electron's magnetic moment is 1000 times that of a nucleon. However, the electron 
polarization is reduced because they are relativistic and degenerate. Furthermore, because of the small v — e cross 
section this polarization may lead to an asymmetry that is smaller than that from nucleon reactions. It is important 
to examine other processes to identify the largest contribution to the asymmetry. 

In this paper, we consider neutrino elastic scattering from slightly polarized neutrons. This may be important (even 
though the neutron polarization is small). Wc do not claim that it is the largest contribution. Instead we focus on 
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elastic neutron scattering for simplicity. The differential cross section is (see for example [14]), 



da/dVL = /^^2 {'^l + ^"^a + {cl - ci,)coiiO + 2P„Ca[(c„ - Ca)cos6'i„ + (c„ + Ca)cos9out]} ■ (1) 

Here G is the Fermi constant, Ei, the neutrino energy (assumed much smaller than the nucleon mass M) and Cy = —1/2, 
Ca = —go/2 with Qa = 1.26. The incident neutrino momentum makes an angle 9in with the polarization direction, 
scatters through an angle 6 and then the outgoing momentum is at an angle Bout with the polarization. 
The polarization of the neutrons P„ depends on the magnetic field B and temperature T, Pn » eB/MT, 

P.»2xl0-lj^l[i^l. (2) 

In the limit = 1 Eq. (1) becomes, 

da/dCl = ao{'i- + PnCosOout), (3) 

with ctq = G'^ E'^ / Att'^ . This simple form provides insight and is good to 10 percent for the asymmetry. [However, we 
use the full result Eq. (1) in the Monte Carlo below.] Equation (3) docs not depend on Oin- Thus the mean free 
path is independent of direction. The asymmetry arises because the outgoing neutrino angular distribution is biased 
towards the polarization direction. 

The total dipole asymmetry A in the neutrino angular distribution I{6) is. 



= y dcose I (9) cose/ j dco&e 1(9). 



(4) 



If the angular distribution from the supernova is proportional to Eq. (3) then A = Pn/3- This asymmetry is related 
to the recoil velocity of the star. 

The gravitational binding energy of a neutron star is of order 100 MeV per nucleon. This is radiated away in 
neutrinos of momentum 100 MeV/n leaving a proto-neutron star of mass about 839 MeV/n. Therefore the recoil 
velocity v of the star is, 

v/c^^A «0.1A (5) 
'839 ^ ^ 

Thus if A is only of order P„ the velocity will be small (around 10~^ of the speed of light c for B near 10^^ G). 

Neutrinos must diffuse through many mean free paths in order to escape the star so they interact repeatedly with the 
polarized neutrons. The crucial question is do these repeated interactions enhance the asymmetry? If the temperature 
distribution at the neutrino sphere is independent of direction than what happens inside may not be important. The 
asymmetry in the neutrino flux will simply arise from the last scattering and be of order P„ . 

However, neutrinos dominate the energy transport. [Note, the effects of convection on A remain to be investigated.] 
Therefore we expect the temperature distribution to be asymmetric because of the asymmetric neutrino flux. This 
could lead to an asymmetry much larger than P„. To investigate multiple interactions, we calculate the asymmetry 
of a "reference configuration" with a very simple Monte Carlo. 

This reference configuration is not meant to be a realistic supernova simulation. Instead it is the simplest system 
with slight nucleon polarization and repeated neutrino interactions. This may allow a simple exploration of the 
physics. We consider a uniform sphere of slightly polarized neutrons. The neutrinos start either at the center or 
uniformly throughout the volume and then interact only via elastic neutron scattering, Eq. (1). We discuss these 
assumptions below. 

The model has two parameters: the polarization of the neutrons P„ and the optical depth of the sphere r/A. This is 
the radius r measured in units of the neutrino mean free path A. For B = 5x 10^^ G and T = 3 MeV the polarization 
is 

P„ » 1 X 10-^ (6) 

Near the center of the star the temperature is larger and the neutrons will become slightly degenerate. This will 
decrease P„ somewhat. Also at high densities, strong interactions may modify P„. It is even possible that there 
is a ferromagnetic phase of dense neutron rich matter [15]. (While this is unlikely at the high temperatures of a 
supernova, there could still be an enhancement in P„.) Perhaps most importantly, we are assuming that the strong 
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B field penetrates the central region. If B is excluded, the average polarization will be lower. Alternatively, there 
could be regions of very high internal fields. For simplicity, we adopt Eq. (6) for P„ and assume that it is uniform 
throughout the sphere. Of course, the final asymmetry is proportional to P„ so it is easy to consider other values. 

We chose the optical depth so that the time scale for neutrinos to diffuse out of the sphere is approximately correct. 
Neutrinos were detected from SN1987A over about 10 seconds [16]. The escape time t is of order the light travel time 
r/c (« 0.1 msec) multiplied by the optical depth r/A. For t to be of order one second requires, 

r/A « 1 X 10^. (7) 

This value is consistent with theoretical simulations [17]. 

We have calculated the overall asymmetry, Eq. (4) using a very simple Monte Carlo code, see Figure 1. To save 
computer time we use a larger P„ = 0.01 than Eq. (6) and smaller r/A than Eq. (7). The results can be scaled to 
the desired values. We find that the asymmetry grows with r/A, 

A^aPn{j) + 0{Pnjf. (8) 

This is our most important result. Note, the second term in Eq. (8) is needed since A saturates for very large r/A. 
The coefficient a arises primarily from the angle averaging in the Monte Carlo. For an r=0 source a w 0.14. For 
uniformly distributed sources a is somewhat smaller a « .057. This refiects the shorter path length for neutrinos 
starting close to the surface. However, A still grows strongly with r/A. 

The linear dependence of Eq. (8) on r/A can be understood with a simple one dimensional biased random walk. At 
each time step the probability to hop left is 0.5 + P„/2 (and right is 0.5 — Pn/2). Because of this bias, the mean value 
of the neutrino's position < a; > is not zero but drifts left with a velocity of order P„. Therefore < x > is proportional 
to t. In contrast, the width of the neutrino distribution grows becaiisc of diffusion but only with t^^"^ . The macroscopic 
asymmetry A depends on the ratio of the mean value to the width. Therefore A grows with t/t^^^ = t^^"^. Finally, 
the time to escape t is proportional to (r/A)^ so A grows linearly with r/A. 

For the parameters of Eqs. (6-7) the reference configuration asymmetry is 

A w 0.014, (9) 

(assuming an r=0 source). This A would imply a recoil of the proto-neutron star of « 400 km/s. This value is 
interesting and much larger than previous estimates. However, it is only a ballpark estimate and we must discuss a 
number of our assumptions. 

First, we only considered neutrino-neutron elastic scattering. For mu and tau neutrinos this is a reasonable first 
approximation to the opacity. However, pair production and annihilation can change the number of neutrinos. 
Multiple interactions will tend to bring neutrinos into thermodynamic equilibrium with the other matter. The 
number of neutrinos then depends on the temperature. Thus, the asymmetric neutrino flux will lead to an asymmetric 
temperature. One side of the star will be warmer than the other side. We expect the asymmetry in the temperature 
to be small near r=0 and grow as one moves out nearer the neutrino sphere. Microscopic simulations to study the 
temperature distribution would be very useful. 

For electron neutrinos one should also consider neutrino capture followed by electron capture. These reactions 
also have asymmetries. For electron capture the asymmetry depends on the electron polarization Pe. which can be 
somewhat higher than for neutrons. However, Pg is multiphed by the small coefficient (c^ — c^)/(c^ + 3c^) « 0.1 [7,9] 
so its contribution may be similar to Eq. (1). 

It is important to discuss antineutrinos. The antineutrino-neutron elastic cross section is related to Eq. (1) by 
switching cosOin with cos0 out, 

da/dCl = ^^2" {^v + + {cl - cl)cosO + 2P„Ca[(c„ + Ca)cos^i„ + (c„ - Ca)cos6out]}- (10) 

This will lead to an opposite sign for the asymmetry. Thus there will be some cancellation between neutrinos and 
antineutrinos. However, this cancellation is not perfect because the opacity is different for electron p and P. Note, 
the physical mechanism is different for u and D. For v, the mean free path is almost independent of direction but the 
outgoing angular distribution is biased. For P, Eq. (3) becomes da/dfl = ao{l + PnCOsOin) so that the mean free path 
depends on direction but the outgoing angular distribiition is almost unbiased. 

Neutrino and antineutrino asymmetries will add in calculating the net lepton number fiux. This could lead to 
significant asymmetries in the trapped lepton fraction and in the chemical composition. If the individual fluxes are 
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say 10 times the net lepton number flux then this may be very sensitive to smafl asymmetries. One percent changes 
in the v and u fluxes will lead to a twenty percent asymmetry in the lepton number flux. This may have an important 
impact on the explosion. For example, on one side of the star the larger neutrino flux may enhance the proton 
fraction while on the other side the antineutrino flux may increase the number of neutrons. One should investigate 
asymmetries in the chemical composition. 

In this paper we examined the asymmetry induced in a supernova from repeated parity violating neutrino-neutron 
elastic scatterings in a strong magnetic field. As a model, we calculated the asymmetry of neutrinos emerging from 
an optically thick sphere of slightly polarized neutrons. We find that the asymmetry grows with the optical thickness 
through which the neutrinos diffuse. This is because the neutrinos interact very many times and leads to parity 
violating effects that are much larger than in previous estimates. 

Future work should examine asymmetries induced in the temperature and chemical composition of the proto-neutron 
star. One should also examine asymmetries in other reactions and the magnetic properties of very dense neutron rich 
matter. We note that parity violation may be sensitive to conditions deep inside the proto-neutron star in addition 
to near the neutrino sphere. In particular, the volume of the core occupied by a strong magnetic field may influence 
the asymmetry. 

This work was supported in part by DOE grant: DE-FG02-87ER40365. 
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FIG. 1. Asymmetry A, Eq. (4), of neutrinos emitted from a neutron sphere of optieal depth r/X. The polarization 
of the neutrons is P„ = 0.01. The sohd curve assumes the neutrinos start at r = while the dotted curve is for 
neutrinos starting uniformly throughout the volume. 
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